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In a recent experiment, the local current distribution of a two-dimensional electron gas in the 
quantum Hall regime was probed by measuring the variation of the conductance due to local gating. 
The main experimental finding was the existence of "hot spots", i.e. regions with high degree of 
sensitivity to local gating, whose density increases as one approaches the quantum Hall transition. 
However, the direct connection between these "hot spots" and regions of high current flow is not 
clear. Here, based on a recent model for the quantum Hall transition consisting of a mixture of 
perfect and quantum links, the relation between the "hot spots" and the current distribution in the 
sample has been investigated. The model reproduces the observed dependence of the number and 
sizes of "hot spots" on the filling factor. It is further demonstrated that these "hot spots" are not 
located in regions where most of the current flows, but rather, in places where the currents flow 
both when injected from the left or from the right. A quantitative measure, the harmonic mean of 
these currents is introduced and correlates very well with the "hot spots" positions. 

PACS numbers: 73.43.-f,73.50.-h 
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I. INTRODUCTION 

The quantum Hall (QH) effect remains a major focus of 
interest^ despite the long time that has passed since its 
discovery. Apparently this is due to the ongoing techno- 
logical progress employing new experimental probes and 
yielding new and sometimes surprising results. A par- 
ticular issue, that has been under debate for quite some 
time, is related to the exact trajectories at which the cur- 
rent flows. Some theories suggest that the current flows 
mainly via edge states along the sample edges^ while 
others, based on the idea of a localization-delocalization 
transition at the centers of Landau levels, predict a distri- 
bution of currents extending throughout the bulk/ How- 
ever, due to the robustness of the conductance quanti- 
zation, the local properties are inaccessible via standard 
transport measurements. While some information on the 
current flow can be derived from static probes^ numer- 
ous attempts have been made to address these questions 
using various scanning imaging techniques^ commonly 
based on local probe of charge and electric potential. Yet, 
while proving successful in describing localized electronic 
states, these methods detect current only indirectly and 
cannot determine unambiguously how the current is par- 
titioned between the edge and bulk channels. 

Recently, a novel experimental approach has been ap- 
plied to probe the local current distribution in a ballistic 
quantum point contact (QPC)i An atomic force micro- 
scope (AFM) tip was placed on top of the 2DEG in which 
a point contact was defined, causing a local depletion of 
electron density beneath it. The underlying assumption 
in this experiment is that this depletion strongly affects 
the conductance through the QPC only if the current 
density under the tip is high. On the other hand, the con- 



ductance should not be modified if there is low current 
density under the tip. Thus, plotting the conductance 
change as a function of tip position results in an imag- 
ing of the electron current density. Indeed, the imaging 
clearly showed the different modes of the electronic wave- 
function being successively occupied as the conductance 
through the QPC increases in quantized steps. Following 
this experiment a theoretical model was devised^ which 
mimics this experiment and yields similar results for the 
distribution of current. 

A similar experimental method has been used more re- 
cently to study the local current distribution in a 2DEG 
in the QH regime^ An AFM tip was placed on top of 
the Hall bar and locally gated the sample, thus chang- 
ing the local potential beneath it. The resistance was 
then measured as a function of tip position and mag- 
netic field. The main finding was that there are "hot 
spots" in the sample, i.e. isolated regions at which the 
conductance is extremely sensitive to the gating poten- 
tial. These domains were interpreted as places where the 
current passes. They are mainly observed in the transi- 
tion between plateaus, and disappear almost completely 
in the quantum Hall regime. 

The fact that these measurement may not directly re- 
flect the total current distribution in the sample can be 
understood through a simple example. Let us assume 
that there is a high barrier separating the left side and 
the right side of the sample. Then a current injected 
from either side (left or right) will be localized on that 
side of the barrier. Nevertheless, the only place where a 
change in the potential may induce a change in the con- 
ductance of the system is at the barrier itself, where no 
current actually flows. To put it in different terms, how 
can the measurement of the conductance, a quantity that 
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obeys specific Onsager relations with respect to reversing 
directions of the current and the magnetic field, yield in- 
formation on the current distribution which, to the best 
of our knowledge, do not obey such symmetry relations? 

Motivated by these experimentally relevant questions 
we employ in this work a recently proposed model for the 
QH transition^ that enables a theoretical modelling of 
the experiment and allows for a detailed analysis of both 
the hot spots and of the spatial distribution of currents in 
the sample. Our main finding, beyond a good qualitative 
agreement with the experimental results, is that the hot 
spots are located at points where current passes both 
when it is injected from the left or from the right, and not 
necessarily at points where the current density is high. 
The number of such symmetry points is enhanced near 
the QH transition, but is small on either side of it, in 
agreement with the experimental results. We propose an 
empirical relation between the current distribution and 
the position of the "hot spots", based on a harmonic 
average of the current distributions when the current is 
injected from the left and from the right. 



II. MODEL 

For the sake of completeness, let us briefly explain our 
model. In strong magnetic fields, electrons with Fermi 
energy tp perform small oscillations around equipoten- 
tial lines. When tp is small, their trajectories are trapped 
inside potential valleys, with weak tunneling occurring 
between adjacent valleys. We associate each such po- 
tential valley with a site in a lattice. Nearest neighbor 
valleys (localized orbits) are connected by links repre- 
senting quantum tunneling between them. As ep in- 
creases and crosses the saddle point energy separating 
two neighboring valleys, the two isolated trajectories co- 
alesce, the electron can freely move from one valley to its 
neighbor and the link connecting them becomes perfect. 
The QH transition occurs when an electron can traverse 
the sample along an equi-potential trajectory (an edge- 
state), which in the model corresponds to percolation of 
perfect links. 

Consequently, the QH problem maps onto a mixture of 
perfect and quantum links on a lattice. Each link carries 
a left going and a right going channel. In accordance with 
the physics at strong magnetic fields there is no scatter- 
ing in the junctions (valleys) and the edge state continues 
propagating uninterrupted according to its chirality (see 
Fig. [IJa)), while the scattering occurs on the link (sad- 
dle point) itself. Each scatterer is characterized by its 
scattering matrix Si, namely a transmission probability 
Ti and phases. The phases are taken as random numbers 
from to 2tt and the transmission amplitude of each link 
is determined by the height of the saddle point barrier 
between the neighboring valleys, taken form a uniform 
distribution U[— V,V]. The transmission is then deter- 
mined locally by the local barrier height e and the Fermi 



energy e F by 



T(e F ) = exp [-a(e - e F ) 



(1) 



for tp < e, where a is some constant, and T(tp) = 1 
for ep > t. The transmission through the system, T, 
is then calculated using the scattering matrix approach 
and the conductance of the system is determined by the 

r, 2 

Landauer formula, G = =f-T. 

Within the scattering model, the current carried by the 
i-th link is given by Ji cx \tp[ 2 — \ip% | 2 , where * s * nc 

incoming (outgoing) wave function from the left (right) 
of the link (see Fig|2J). Note that due to the unitarity of 
the S-matrix, the current is a locally left-right symmetric 
quantity, and is thus a property of the entire link. 




FIG. 1: Mapping of the problem onto a quantum percola- 
tion on a lattice: each link carries two counter propagating 
edge modes (a). A non zero transmission (thin lines) allows 
electrons to tunnel between adjacent sites (potential valleys) . 
When the transmission is unity (bold lines in (a)), these two 
valleys merge, and an edge state can freely propagate from 
one to another. A percolation of these perfect transmission 
links (b) correspond to an edge state propagating through the 
system and a quantized conductance. 
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FIG. 2: Within the scattering model, each link is character- 
ized by a scattering matrix Si, connecting between the in- 
coming and outgoing wave functions from left to right (see 
text). 



As the Fermi energy increases, the conductance rises 
from zero to unity at the percolation threshold (Figs.^b) 
and[3J). In Ref. 14 it was demonstrated that the phase 
transition described by the above model exhibits a criti- 
cal exponent v ~ 2.4, in agreement with numerical sim- 
ulations for other models describing the QH transition. 
In this paper we simulate the experiment of Ref. |9| using 
the above model. 
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FIG. 3: Conductance as a function of probability p. The ar- 
rows point on the concentrations at which the current distribution 
(Fig-EJ an d the hot spots (Fig. 1101 are plotted. 



III. RESULTS 



Current Distribution 



tween the distribution of left- and right-coming currents. 
A similar study, based on a tight-binding description, 16 
have resulted with similar spatial current distributions, 
and has emphasized the role of the current chirality in 
the quantization of the Hall conductance and the vanish- 
ing of the longitudinal resistance far from the transition. 
The left-right a-symmetry in the presence of magnetic 
fields was also pointed out in Ref. 
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In Fig. |31 the conductance (in units of 2e 2 /h) of a 
20 x 20 size system^ is plotted as a function of prob- 
ability of a link to be perfect which corresponds to Fermi 
energy, p = ^ + l). As can be seen, our model re- 
produces the conductance fluctuations close to the perco- 
lation transition^ , which occurs (for this realization) at 
p = 0.505, close to the bulk classical percolation critical 
point. 

The current distributions obtained for a specific real- 
ization of disorder for different values of p are depicted 
in Fig.^lfor two cases, when the current is injected from 
the left (left column) and from the right (right column), 
bright colors correspond to a high current. Each row 
corresponds a concentration p as denoted by arrows in 
Fig. |3 As can be seen, for low concentrations (i.e. on 
the insulating side of the transition, Fig. CHa)-(b)) the 
currents flow in the system along some closed trajectory, 
returning back to the lead they came from due to the 
low transmission. In this realization the potential bar- 
rier separating left and right is located closer to the left 
side of the sample, thus an electron injected from the 
left is reflected almost immediately, while an electron in- 
jected from the right meanders through a larger part of 
the sample before being reflected. 

Close to the transition the amount of current that 
passes through the system from one side to the other is 
roughly equal to the amount of back-scattered current, 
since the transmission of the system is close to T rj 0.5 
(Fig. E3c)-(d)). Here we find that the correlation be- 
tween left- and right-originating current distributions is 
higher, and that the spatial distributions are broad, in 
accordance with percolation theory. Finally, for large 
concentrations (Fig. EJe)) the transmission of the sys- 
tem is perfect, and thus the current passes between the 
leads following some trajectory, which corresponds to a 
percolating path of perfect links. As in the case of low 
concentration, the chirality causes strong separation be- 
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FIG. 4: Spatial current distribution when current is injected from 
the left (left column) or from the right (right column), for the con- 
centrations depicted by arrows in Fig-El I n t ne numeric calculation 
the current is injected from the upper-left-most or lower-right-most 
link, corresponding to left- or right-injected current. As seen, while 
far from the transition left- and right-coming currents are spatially 
separated, close to the transition the percolative nature of the sys- 
tem causes a wide spatial current distribution. 
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B. Hot Spots 

In order to simulate the "hot spots" experiment^ we 
mimic the local gating at a certain site by adding addi- 
tional energy to the potential barriers surrounding that 
site, up to a distance of several lattice spacings. The 
conductance of the original lattice is then compared with 
that of the perturbed lattice for different values of the 
Fermi energy (corresponding to the experimental change 
in the magnetic field, or filling factor). 

In Fig. the spatial distribution of the change in 
the conductance is plotted for different values of p, cor- 
responding to the points denoted 6, c, d, e in Fig. [3] 
Brighter points correspond to the " hot spots" of the cur- 
rent, for which there is a sizable change in the conduc- 
tance as that point in the lattice is gated by the AFM 
tip. The numerical data are normalized to the largest 
conductance change. 



(a) p=0.39 (b) p=0.42 




(c) p=0.465 (d) p=0.5025 




FIG. 5: Spatial mapping of the "hot spots" in the conductance, 
for probabilities corresponding to the arrows (b)-(e) in Fig.|3 Far 
below the QH transition (a) there are very few hot spots. Their 
number and intensity increases as one approaches the transition 
(b), is maximized at the transition region (c) and again diminishes 
above the transition (d). 

As deduced from Fig. both the amount and inten- 
sity (i.e. how considerable the influence of depleting sites 
from electrons on the conductance is) of the hot spots 
increase as one approaches the QH transition, in agree- 
ment with the experimental result. This may be seen 
more clearly by calculating the absolute value of the con- 
ductance change 5a, averaged over all the lattice sites 
and over disorder. 5a is plotted in Fig. for 100 disor- 
der realizations. It is found that the maximal change in 
the conductance corresponds to the percolation thresh- 



old, determined by the point at which the conductance 
is length-independent, denoted by an arrow in the inset 
of Fig. El 
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FIG. 6: Conductance (stars) and absolute value of the change 
in the conductance (triangles) as a function of concentration p, 
averaged over all lattice sites and over 100 realizations. The arrow 
indicates the point of percolation, denned as the point at which 
the disorder-average conductance is independent of system length. 
The inset shows the conductance for system lengths L = 10, 15, 20. 



C. Effect of tip parameters 

Let us turn our attention to the role of the AFM tip. 
Experimentally, the exact effect of the AFM tip on the 
sample is unknown. However, it is reasonable to assume 
that the tip induces an increase in the potential energy 
in the area underneath it, which affects electrons up to a 
length l e ff away from it. 

In order to address this point theoretically, we note 
that the tip has two main tunable parameters, namely 
the potential difference (voltage) between the tip and the 
sample, Etip, and the distance between the tip and the 
sample. Changes in these experimental parameters af- 
fect two different aspects of the model: (a) a change in 
the offset of the energy in the links underneath the tip, 
namely the local potential energy change induced by the 
tip, and (b) a change in the effective length, l e g, over 
which electrons feel the tip. While experimentally these 
two parameters are both affected, to some degree, by 
the tip height and voltage, theoretically one can study 
the change in each parameter separately. To simulate 
these effects, we repeat the above calculation, with a tip- 
induced exponential-enveloped change in the local poten- 
tial barrier on the links, AE — Eti p e~ d / lc!! , where d is 
the distance between the link and the position of the tip. 
In what follows we explore how changing either Ef lp or 
l e g affects the conductance change. 

Changing E t i p does not have a significant effect on the 
spatial distribution of the hot spots, but only on their 
strength, namely the conductance change induced by the 
tip. In Fig. [3 we plot the average change in conductance 
5a (averaged over the entire sample) as a function of -Etip) 
varied from E tlp = to E t - lp = V, that is of the order 
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of the band width. The calculation is performed for con- 
centration p — 0.42125, where the number of hot spots is 
quite large, and i e fl = 1 ( m units of lattice spacing) . As 
seen, 5a increases monotonically with -Etip- In the inset 
of Fig. we plot the conductance change 5a when the 
tip is over the strongest "hot spot" in the sample, and 
again, a monotonic increase in 5a is observed. 
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FIG. 7: Average change in conductance, <5<r, as a function of STM 
tip voltage -Etip, averaged over the entire sample. Inset : the same 
for the strongest "hot spot" in the sample. 



Next we examine the effect of changing Z c ff . One may 
naively guess that increasing should result in an in- 
crease in 5a. However, due to quantum interference this 
may not always be the case, especially far from the perco- 
lation transition, where the conductance change is rather 
small. In Fig. |H1 We plot the spatial map of hot-spots 
for i^tip = V and p — 0.30875 (which is close to the 
percolation transition for this sample, the sample hav- 
ing a transmission T = 0.648 for this concentration), for 
Z c ff = 1, 2, .., 8 (in units of lattice constant). We find that 
although the number of hot-spots (indicated by bright 
colors) increases, their locations are changed, due to in- 
terference effects. 

To make this more qualitative, in Fig. we plot the 
average change in the conductance 5a, averaged over the 
sample and over 100 realizations of disorder, for concen- 
tration p = 0.5 (i.e. with Fermi energy at the center of 
the band), as a function of l e s . We find that indeed, on 
average the conductance change exhibits a monotonic in- 
crease. However, for a given realization and a given tip 
position, changing Z c ff (which corresponds to a change 
in the distance between tip and sample) may result in 
fluctuations in 5a, as plotted in the inset of Fig. E3 




FIG. 8: Spatial map of hot-spots for different values of the STM tip effective length, i e g = 1,2, ..,8. 



D. Relation Between Hot Spots and Current carrying electronic states in the bulk ? As stated above, 

Distribution although one is inclined to give a positive answer, the 

different symmetry of the hot-spots and the current dis- 

Next, we ask the question: are the "hot spots" ob- 
served in experiment located at the extended, current 
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FIG. 9: Average conductance change, <5cr, averaged over the sample 
and over 100 realizations of disorder as a function of Z cf j , showing a 
monotonic increase. Inset : the local change 8<x at a given position 
and realization of disorder as a function of Z c ff . One sees that <5<r 
may fluctuate due to quantum interference. 



tribution points that they cannot be identical. It is clear, 
however, that they are correlated, in a way we discuss be- 
low. 

Let us examine the correlation between the hot-spots 
and current distribution 17 . Since the strength and loca- 
tion of the hot links is independent of the direction of 
the current injection, it is clear that the hot links are not 
located at the points where most of the current passes, as 
these depend sensitively on direction. Rather, it is found 
that the hot links are located at points which hold ap- 
preciable currents for both directions of current injection. 
To demonstrate this we plot on the left column of Fig.llUI 
the spatial image of a harmonic average of the local cur- 
rents from the two directions of current injection. Bright 
spots thus correspond to links in which local current is 
significant in both directions of current injection. On the 
right column we plot the spatial image of the "hot links" , 
that is, the conductance change due to setting the trans- 
mission on each link to zero. Both columns are shown for 
the concentrations corresponding to the arrows of Fig |3] 
One clearly sees the correlation between the two images. 

In order to make this correlation more quantitative 
the correlation function C (p) is defined to be the square 
of difference between the current and the conductance 
change, normalized and averaged over the entire sample. 
That is, let ji be the normalized current in the «-th link 
(either when the current is injected from the left, from 
the right or a average of the two), and 5a \ be the nor- 
malized change in the conductance when the tip is placed 
over the i-ih link (note that both ji and Sai depend 
on the concentration p), then C(p) — jj J2i \h ~ ^ a i\ 2 - 
The smaller C(p) is, the higher is the correlation be- 
tween these distributions. In Fig. El C(p) is plotted for 
different concentrations p, when correlating between the 
hot links and the harmonic averaged current distribution 
(triangles), the hot-links and the left-coming (stars) and 
right-coming (squares) current distribution and between 
the hot-links and a random distribution (diamonds), that 
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FIG. 10: Left column: a spatial image of the hot links, for the con- 
centrations depicted by arrows in Fig. [3] bright links correspond to 
links with strong effect on the conductance. Right column: spatial 
image of the harmonic mean of local currents obtained from two 
different directions of voltage drop, for the same concentrations. 
Bright links correspond to links in which current is considerable 
for both directions of voltage drop. A clear correlation between 
the left and right columns is visible. 



serves as a reference scale. As seen, the hot-links are well 
correlated with the harmonic-averaged current distribu- 
tion, indicated by the low values of C(p). The correlation 
between the hot-links and currents flowing from the left 
or from the right is much worse, and actually resembles 
the correlation with a completely random distribution for 
large concentrations. 
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FIG. 11: Correlation function C(p) (see text), as a function of 
concentration. The correlation is between the hot-links and the 
harmonic-averaged current distribution (triangles), the hot-links 
and the left-coming (stars) and right-coming (squares) current dis- 
tribution and between the hot-links and a random current distri- 
bution (diamonds). The lower the value of C(p), the higher the 
correlation. The low value of the correlation function between the 
hot spots and the harmonic-average, relative to the other correla- 
tions, demonstrate that the hot spots do not directly reflect the 
total spatial current distribution, but only the transport current, 
estimated by the harmonic average. 

IV. SUMMARY 

In this work a recent model for the QH transition 14 
have been used to shed light on an experiment in which 
a local AFM-tip induced local gating has been employed 



to study the change in the resistance as a function of the 
tip location. The present theoretical study demonstrates 
the existence of " hot spots" , regions with higher sensitiv- 
ity to local gating. Our model, consisting of a coherent 
mixture of perfect and quantum links, qualitatively re- 
produces the experimental result. It was demonstrated 
that these "hot spots" do not lie in the current-carrying 
paths, but rather on areas in the sample where current 
flows both when it is injected from the left or from the 
right, that is on the left-right symmetric parts of the 
current carrying paths. Note that since the geometry 
employed in this study corresponds to a two-terminal ge- 
ometry, the longitudinal conductance is trivially related 
to the Hall conductance, and thus one expect a similar 
behavior of the "hot spots" in the Hall conductance. 

We conclude by noting that in order to verify our find- 
ing experimentally, one should imply a non-destructive 
local probing of the QH sample (that is, probing that does 
not affect the conductance), e.g. local current-induced 
magnetic field sensing, in addition to the above men- 
tioned local AFM-tip gating, on the same sample. Such 
an experimental system may turn out to be useful for 
measuring local currents in other systems, such as disor- 
dered superconducting thin films, or Hall samples in the 
fractional QH regime. 

We thank T. Ihn for valuable discussions. This re- 
search is partially supported by a grant from the Israeli 
Science Foundation (ISF). 
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